Making Quantum Spin-Hall Effect Robust via Magnetic Manipulation 
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The quantum spin Hall (QSH) effect is known to be unstable to perturbations violating time- 
reversal symmetry. We show that creating a narrow ferromagnetic (FM) region near the edge of a 
QSH sample can push one of the counterpropagating edge states to the inner boundary of the FM 
region, and leave the other at the outer boundary, without changing their spin polarizations and 
propagation directions. Since the two edge states are spatially separated into different "lanes" , the 
QSH effect becomes robust against symmetry-breaking perturbations. 
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The quantum spin Hall (QSH) effect, a quantum state 
of matter, has attracted much attention in recent years, 
because of its fundamental interest and potential ap- 
plications in spintronic devices. The QSH effect was 
first predicted theoretically by Kane and Mele [H and 
by Bernevig and Zhang Q, in independent works, and 
soon after observed experimentally in HgTe quantum 
wells. [1, 01 A key ingredient to the QSH effect is a 
strong intrinsic spin-orbit coupling, which acts as spin- 
dependent magnetic fluxes coupled to the electron mo- 
mentum. In the ideal case, where electron spin is con- 
served, the two spin sectors of a QSH system behave 
like two independent quantum Hall (QH) systems with- 
out Landau levels. [5] They contribute opposite quantized 
Hall conductivities, when the electron Fermi level is in- 
side the bulk band gap, so that the total Hall conductiv- 
ity vanishes but the spin Hall conductivity is quantized. 
On a sample edge, two counterpropagating gapless edge 
modes with opposite spin polarizations exist in the bulk 
band gap, which can transport spin currents without dis- 
sipation of energy. 

When the spin conservation is destroyed, e.g., by the 
Rashba-like spin-orbit coupling, the spin Hall conductiv- 
ity deviates from the quantized value. 0] However, the 
edge transport can remain to be dissipationless, pro- 
vided that the time-reversal (TR) symmetry is present 
and the bulk band gap is not closed. In this case, a QSH 
system can no longer be divided into two QH systems, 
and the existence of the gapless edge states has been at- 
tributed to the nontrivial topological properties of bulk 
energy bands. The nontrivial bulk band topology of the 
QSH systems is usually described by the Z2 index Q or 
the spin Chern numbers. These topological invari- 



ants reveal the fundamental distinction between a QSH 
insulator and an ordinary band insulator. 

Extension of the idea of the QSH effect to higher di- 
mension has led to the proposals of three-dimensional 
(3D) TIs. Research of the 3D TIs has been fruit- 



ful in recent years, both theoretically and experimentally. 
A 3D TI has a bulk band gap and gapless surface states 
on the sample boundary. A simple 3D TI consists of lay- 



ers of 2D QSH systems. However, such a 3D TI is not 
stable in the presence of disorder due to interlayer scat- 
tering, and so called the weak TI. A strong TI is topolog- 
ically nontrivial in any direction, and the surface states 
are characterized by an odd number of Dirac cones. The 
metallic surface states of the 3D TIs provide a unique 
platform for realizing some exotic physical phenomena, 
such as Majorana fermions [17[ and topological magneto- 
electric effect. 18, The existence of topological surface 
states in the 3D TIs has been experimentally confirmed 
in Bii-xSba;, Bi2Te3, and Bi2Se3 materials, '2d-25] which 
evokes a great surge of research interest in this field. 

While the TR symmetry was often considered to be a 
prerequisite to the QSH effect, its role is two-sided. In 
a TR invariant QSH system, the two oppositely moving 
edge states at the Fermi energy are connected to each 
other under TR, and so have opposite spin orientations. 
Elastic backscattering from nonmagnetic random poten- 
tial is forbidden. On the other hand, the two opposite 
movers have identical spatial probability distributions. 
Turning on small TR-symmetry-breaking perturbations 
immediately couples the two edge states, giving rise to 
backscattering. This makes the QSH effect fragile in real- 
istic environments, where perturbations violating the TR 
symmetry are usually unavoidable. Experimentally, two- 
terminal conductance close to the predicted value 2e^//i 
was observed only for small QSH samples with dimen- 
sions of about (1 X l)^rn? [3], in contrast to the tra- 
ditional QH effect, where the Hall conductivity can be 
precisely quantized on macroscopic scales. So far, QSH 
effect as robust as the QH effect has been elusive. It was 
found recently that the nontrivial bulk band topology of 
the QSH systems remains intact, even when the TR sym- 
metry is broken, '26^ implying that the instability of the 
QSH effect is solely due to properties of the edge states. 

In this Letter, we show that the QSH effect can be 
stablized by designed removal of the TR symmetry of 
the edge states. Creating a narrow ferromagnetic (FM) 
region near one edge of a QSH sample can push one of 
the counterpropagating edge states to the inner bound- 
ary of the FM region, and leave the other on the outer 
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boundary. The moving directions and spin orientations 
of the individual edge states are unchanged, so that the 
QSH effect persists. Importantly, the edge states are spa- 
tially separated, so that the QSH effect becomes robust 
against general perturbations without fictitious symme- 
try constraints. Our work paves a road to realize robust 
QSH effect via magnetic manipulation. 

We start from the effective Hamiltonian for a HgTe 
quantum well with a uniform FM exchange field, given 
by 

H = v^[T^Ka^ + ky&y) + (Mq - Bk'^)a,, + Hi , (I) 

where r stand for Pauli matrices for two spin states, 
and (T for the electron and hole bands. An exchange 
field can be created in the Hg Te layer by doping of mag- 
netic atoms, such as Mn, [27| or through magnetic prox- 
imity effect with a FM insulator. The exchange field 
can be described by (27| Hi = goTzO'z + gA, where 
go = ^(G'h-Ge) and.gi = ^{G^+G^) with 26^ (2Gh) as 
the exchange splitting of the electron (hole) bands. For 
convenience, we set the reduced Planck constant h to be 
unity, and take the momentum to be dimensionless by 
properly choosing the units for parameters of the model, 
namely, v^, Mq, B, go, and gi. 

It is worth mentioning that the model Hamiltonian Eq. 
(f ) can describe not only the HgTe quantum wells, but 
also a class of QSH systems. [28| For example, setting 
gi = and making a unitary transformation "H — U'' HU 
with U = + + (1 - f^)ay]e-''T'^', we obtain V. = 
Vp{kyax — kxCTy) + (Mq — Bfc^) t^ctz +goO'z- One can find 
that this Hamiltonian is identical to that of a thin film of 
3D TI Bi2Se3 with an exchange field, 11, 2^ for which 
<T stand for spin up and down, and r for bonding and 
antibonding of the surface states on two surfaces. 

Since is a conserved quantity, one can easily diag- 
onalize Eq. (1), and obtain two conduction bands and 
two valence bands. Under the condition of \go\ > \gi\, 
a nonzero middle band gap exists except at go = iAfo, 
where the conduction and valence bands touch. As has 
been discussed in Ref. [23|, for HgTe quantum wells 
doped with Mn, Ge and Gh have opposite signs, so that 
the condition above is always satisfied. The spin Chern 
numbers for T2 = ±1 can be derived to be 

G± = ±i[sgn(i?) + sgn(Afo ± 50)] • (2) 

At 50 = 0, G± = ±sgn(i?) if BMo > 0, corresponding to 
a QSH phase, and C± = otherwise, corresponding to an 
ordinary insulator. We focus on systems with BAIo > 0, 
which are in the QSH phase at go = 0, and for the sake 
of definiteness, we will confine ourselves to the parameter 
region of S > and Mq > 0, in which C± ~ (1, —1) at 
go = 0. (All the conclusions reached in this work do also 
apply to the case of B < and Mq < 0.) Notably, with 
increasing go to go = Mo, C± undergo a transition from 
(1,-1) to (1,0).' 




FIG. 1: Schematic of the QSH sample with a long strip geom- 
etry. The profile of the y-dependent exchange field is shown 
by the thick (red) line. 



Next we consider a QSH sample with a long strip ge- 
ometry, as shown in Fig. 1. We will set gi = for now, 
which does not affect topological properties of the sys- 
tem in the case of \go\ < \gi\. The effect of finite gi will 
be discussed later. The exchange field go is taken to be 
nonzero in narrow regions of width I near the two edges 
of the strip, namely, region I of < y < Z and region 
II oi {L — I) < y < L, and vanishing in region HI of 
I < y < (L — l). The system as a whole has a bulk energy 
gap around energy 0, for go ^ ±Mo. The edge states in 
the bulk energy gap can be solved analytically by replac- 
ing ky with ~idy in the system Hamiltonian. For = 1, 
an edge mode with energy E+{kx) — found on 

one side of the strip, whose wavefunction is given by 



V+{kx,y) = |1, -l)(t>+ikx,y) . 



(3) 



Here, ket \Tz,ax) with = ±1 and = ±1 is used to 
represent the common eigenstate of and ax- The spa- 
tial wavefunction (l>+{kx,y) = Cie-J'/^i^s) - e-^'/^^^f'] for 
y <l, and (t)+{kx,y) = L»ie-(?^-''/«i(o) - D2e-'^y-^y^^^^°^ 
for y > I, with G, Di and D2 to be determined from the 
conditions of continuity and normalization of (j)+{kx,y)- 
The two characteristic length functions are defined as 



Ci,2(e) 



2B 



Vp ± - 4B(Mo - Bkl + e) 



(4) 



For Tz = —1, we find another edge mode with energy 
E-{kx) = vpkx and wavefunction 



(p-{kx,y) = I - l,-l)<p-{kx,y) , 



(5) 



where (t>-{kx,y) = E[e-y/^^^-3) - e-y/^^^-^)] for y < 
I, and 0-(fcx,y) = Fie-(s'-')/«i(o) _ Fse-^^'-'^/^^t") for 
y > I. Owing to the two- fold rotation symmetry, the 
edge modes on the other side of the strip have dispersion 
relations E±{kx) — -iiv^kx- Their wavefunctions can be 
obtained through replacements y L — y, Ux ^ —^x 
and kx — > —kx in Eqs. ^ and (HJ. The energy spectrum 
of edge states appears to have mirror symmetry between 
kx and —kx- However, this is an accidental event, and 
the system does not possess mirror symmetry if go ^ 0. 

In Eqs. ^ and JS]), the r and a parts of wavefunc- 
tions do not change with varying go- In Fig. 2(a), the 
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FIG. 2: (a) \(j>+(kj:,y)f and |(^_ (fco;, at fc^r = as func- 
tions of y for different go, where Vp = 1, B — 5, Mo = 0.2, and 
I = 100 are taken. Here, \(l>+{kx,y)\^ is plotted only for two 
values of go, as its change with go is relatively small. Spatial 
distributions of spin Chern numbers C± and edge states are 
shown for go < Mo (b) and go > Mo (c). 
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FIG. 3: Upper panel: calculated energy spectra for four sets 
of go and 7. Lower panel: magnitude of the energy gap of 
edge states as a function of go for three values of 7. The 
other parameters are taken to he Vf = 1, B = 5, Mo = 0.2, 
L = 1000, and I = 200. 



modulus squared of spatial wavefunctions 0_|_(O,y) and 
0„(O,y) are plotted as functions of y for several values 
of go- The Fermi velocity Vjr is taken to be the unit of 
energy. At 50 = 0, we have \(j)+{0,y)\'^ = |0_(O,y)p, and 
both lines coincide with each other, which is required by 
the TR symmetry at go = 0, as mentioned above. With 
increasing go, the peak of |0+(O,y)p becomes sharper 
and closer to y = 0. On the contrary, the shape of 
|(/)_(0,?/)p widens with increasing go, and spreads across 
region I at go = Mo — 0.2. With further increasing 
go, |^_(0,?/)p becomes localized near y ^ I ^ 100, i.e., 
the inner boundary of region I. This transition is di- 
rectly caused by the sign change of £,2{—g) in 0-(O,y) 
at 5o = Mo, from positive to negative. The other charac- 
teristic lengths in (/)_(- (0,j/) and (/)_(0,y) remain positive 
during this process. 

The evolution of the edge states with varying go is 
further illustrated in Figs. 2(b) and 2(c), and can be un- 
derstood in terms of calculated spin Chern numbers. For 
go < Mo, the spin Chern numbers C± in the three regions 
take the same value (1,-1). This indicates that the three 
regions are topologically equivalent, and can be regarded 
as a QSH system as a whole. As a result, the edge states 
for both up spin (r^ = 1) and down spin (tz = —1) ap- 
pear near the sample boundaries y = and y = L, as 
shown in Fig. 2(b). For go > Mo, the situation is quite 
different, because C± undergo a transition at go — Mo 
from (1,-1) to (1,0) in regions I and II. For the spin- up 
electrons, the three regions have the same Chern number 
C+ = 1, and as a whole are equivalent to a QH system. 
Therefore, the spin-up edge states remain localized near 
y = and y = L. For the spin-down electrons, region III 
with C- = — 1 is a QH system, sandwiched between two 



insulators in regions I and II, where C_ = 0. The spin- 
down edge states thus shift to their interfaces, namely, 
y — I and y — L — I, as shown in Fig. 2(c). 

Comparing Fig. 2(c) for go > Mo with Fig. 2(b) for 
go < Mo, one finds that both systems have very similar 
edge states: the lower edge contains a left mover with up 
spin and a right mover with down spin, and conversely 
for the upper edge, so as to exhibit the same QSH ef- 
fect. An important difference is that for go > Mo, the 
counterflows of electrons at the lower (upper) edge are 
spatially separated into two different "lanes" located at 
y ^ (y ~ L) and y = 1 {y — L — I), which provides an 
essential protection. [30| As a result, the QSH state for 
go > Mo should be more stable than that for go < Mo- 
To show this point, we introduce a specific perturbation 
violating the TR symmetry. Since the two edge states 
given by Eqs. (3) and (5) have opposite r^, we include 
an in-plane Zeeman field throughout the sample with po- 
tential 

V^lfy, (6) 

which represents a "maximal" coupling between the two 
edge states near the degenerate point. A tight-binding 
model on a square lattice is constructed for the strip ge- 
ometry, which recovers the form of Eq. (1) in the con- 
tinuum limit. The energy spectrum is calculated numer- 
ically from the tight-binding model. 

The calculated energy spectra for four sets of go and 
7 are plotted in Fig. 3. Since the dispersion relation for 
edge states is given by E±{kx) = ^Vj^kx on one side of 
the strip and E±{kx) = iv^kx on the other side, the 
energy spectra of edge states on the two sides collapse, 
and only two branches of edge states are seen in each 
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FIG. 4; Calculated energy spectra with 7 = 0.01 and gi — 
0.3(70 for go = 0.1 (a) and go — 0.4 (b), and spatial proba- 
bility distributions of the two edge states near the lower edge 
at energy E = —0.01 for go = 0.1 and 0.4 (c). The other 
parameters are the same as in Fig. 3. Here, spin up and spin 
down refer to the edge states with positive and negative spin 
polarizations (fz), as Tz is not conserved. 



breaks the particle-hole symmetry. It is found that the 
edge states are gapped for 50 = 0-1 (< -^o) and become 
gapless for go = 0.4 (> Mq), respectively, shown in Figs. 
4(a) and 4(b). The spatial probability distributions of 
the two edge states at E — —0.01 on one side of the strip 
are plotted in Fig. 4(c). While the spin- up polarized edge 
state stays near y = 0, the spin-down polarized edge state 
moves from near y = to y = l = 200 with tuning go 
from 0.1 to 0.4, exhibiting the same transition as in Fig. 
2(a). This result indicates that the presence of nonzero 
gi does not change the physical picture qualitatively, and 
the robust QSH effect can be achieved when the exchange 
field is sufficiently strong. 
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figure. For go — j — 0, the edge states are gapless, as 
the system is TR invariant, and is conserved as well. 
Upon turning on 7 = 0.01, an energy gap appears in the 
edge state spectrum even at go = 0, signaling the onset of 
backscattering. With increasing go to 0.1, the energy gap 
decreases but remains finite. With further increasing go 
to 0.4, greater than Mq = 0.2, the energy gap vanishes, 
an indication of quenching of backscattering. The mag- 
nitude of the energy gap of the edge states is plotted as 
a function of go in the lower panel of Fig. 3, for three 
different values of 7. While it increases almost propor- 
tionally to 7, the energy gap decreases with increasing 
go- When go slightly exceeds the critical point Mq = 0.2, 
the energy gap essentially vanishes for all the three values 
of 7. Furthermore, we find that for go < Mo, the spin 
polarizations of the edge states (fz) deviate from ±1 ap- 
preciably near the energy gap, as nonzero 7 destroys the 
conservation of f^, but nearly full spin polarizations are 
restored for go > Mo- These numerical results confirm 
our argument that the QSH effect in the go > Mq region 
is robust against general perturbations without symme- 
try constraints. 

Finally, we investigate the effect of nonzero gi in the 
Hamiltonian on the QSH effect. Nonzero gi together 
with go is present in FM regions I and II of Fig. 1. The 
calculated energy spectra are plotted in Fig. 4 for go = 
0.1 and 0.4, with 171 — 0.3.go set in proportion to go. The 
energy spectra of edge states on the two sides of the strip 
are no longer degenerate in the presence of nonzero gi, 
the positions and spin orientations of the split edge states 
being labeled in Fig. 4(b). Inclusion of nonzero gi also 
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